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<£ ■ Abstract 

m 

We study in some detail the next-to-minimal super symmetric standard model 
with gauge mediation of supersymmetry breaking. We find that it is feasible 
qq \ to spontaneously generate values of the Higgs mass parameters \i and B^ 

consistent with radiative electroweak symmetry breaking. The model has a 



phenomenologically viable particle spectrum. Messenger sneutrinos with mass 



C i n the range 6 to 25 TeV can serve as cold dark matter. It is also possible to 



evade the cosmological domain wall problem in this scenario. 
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I. INTRODUCTION 

The Standard Model (SM) fulfills electroweak symmetry breaking (EWSB) by introduc- 
ing a scalar Higgs particle. The mass-squared of this particle suffers from quadratically 
divergent radiative corrections, and is therefore sensitive to physics effects in the ultraviolet 
regime, possibly at the grand unification scale (Mqut) or the Planck scale (Mpi). Stabilizing 
the electroweak scale against radiative corrections is the primary motivation for exploring 
theories beyond the SM, among which weak-scale supersymmetry (SUSY) has been the 
leading candidate 0. 

In the minimal supersymmetric extension of the standard model (MSSM), the Higgs 
scalar potential has the form 

VD(fi 2 + m\ u ) \E U \ 2 + (/i 2 + m% d ) \H d \ 2 - {B, t H u H d + h.c.) + ... (1) 

where H u and H d are two Higgs doublets coupled to the up and down type fermions respec- 
tively, rriH u d and B^ are soft SUSY breaking parameters, and fi is a mass parameter in the 
superpotential 

W» D fiH u H d . (2) 

To achieve correct electroweak symmetry breaking without drastic fine-tuning, all these mass 
parameters should have the size of 0{Mz). The soft parameters can be made to satisfy this 
requirement. The \x parameter on the other hand, should naively be of the order of the 
fundamental scale of the theory like Mpi or Mqut since it is a dimensionful parameter of 
the superpotential. This apparent mismatch is termed as the SUSY yU-problem. 

There are many approaches to resolving the /z-problem |2|||, the simplest one being the 
introduction of an extra singlet Higgs field N. The Higgs superpotential is then written as 

W N D \ N NH u H d - ^N 3 . (3) 

This is called the next-to-minimal supersymmetric standard model (NMSSM) ||. The 
bilinear term iiH u H d is absent by invoking a discrete Z3 symmetry under which every chiral 
superfield $ transforms as $ — > e 27n//3 $. When N acquires a vacuum expectation value 
(VEV) (N) ~ O(Mz), the /i-parameter is effectively generated as 

fx = \ N (N) ~ 0(M Z ). (4) 

The inclusion of the iV 3 term in Eq. (§) is important since without it, the superpotential 
would have a problematic Peccei-Quinn symmetry. The spontaneous breaking of the Z3 
symmetry by (N) leads to the well-known cosmological domain wall problem ||. Also, one 
should be aware of the potential problem of destabilizing the electroweak scale from tadpoles 
associated with a gauge singlet [|J. 

In gauge-mediated supersymmetry breaking models [§-Qj, the SUSY breaking in the se- 
cluded sector is usually parameterized by a SM singlet chiral superfield S with non-vanishing 
VEVs for the lowest and highest components (S) + 6 2 Fs- If the superpotential contains the 
term 



W s D \ s SH u H d . (5) 

then /i = Xs(S) and B^ = A5F5, which implies 

Bp ~ (iA, (6) 

where A = F S /(S) ~ 0(100 TeV). Therefore \i and ^B^ cannot both be 0(M Z ). This 
is known as the /z-problem in GMSB. Some solutions have been proposed in PJH|. In 
particular, Ref. |1(J studied the GMSB paradigm with the NMSSM, but the conclusion is 



negative if one does not introduce extra vector-like quarks in the theory. 

In this paper we perform a detailed analysis of a gauge-mediated model with an extra 
singlet Higgs field in a most general superpotential including all couplings between the 
singlets, messengers and Higgs fields, respecting the Z3 symmetry. Without constructing 
a SUSY breaking model explicitly, we consider it as a phenomenological model and study 
its consistency with theoretical and experimental constraints. In Section II we describe our 
general framework for the gauge-mediated SUSY breaking model with an extra singlet. In 
Section III we present a numerical study of this model, concentrating on the /i-problem and 
phenomenologically viable mass spectrum of SUSY particles (sparticles) with emphasis on 
the next-to-the-lightest supersymmetric particle (NLSP) and Higgs bosons. In contrast to 
the conclusion of Ref. |L(J, we find that suitable /i and B^ parameters can be generated 
in this minimal scenario with all phenomenological constraints satisfied. In Section IV we 
show that messenger sneutrinos can serve as a cold dark matter candidate and that they 
are sufficiently massive to evade the current experimental bound from direct searches. In 
Section V we demonstrate that the domain wall problem may find a solution due to fast 
decay of the walls induced by higher dimensional operators, yet consistent with EWSB. 
We conclude in Section VI. Some technical details of the model — the renormalization group 
equations (RGEs) and Higgs boson mass matrices — are provided in two appendices. 

II. THE MINIMAL GMSB WITH AN EXTRA SINGLET 

In our phenomenological approach to the GMSB, we parameterize SUSY breaking in the 
secluded sector by a SM gauge singlet spurion S = (S) + 2 F S . The SUSY breaking is then 
transmitted to the observable sector via vector-like messengers $ and $ which couple to S 
according to 5*$$. In the simplest case, $ and $ transform as a single flavor of 5 + 5 of 
SU(5). 

In gauge-mediated models not involving any direct messenger-matter interactions, it is 
not possible to generate fi and B^ consistent with EWSB. In fact, even if the spurion 5* 
couples directly to the Higgs superfields as in Eq. (|[), one can easily show that fi = Xs{S) 
and B^ = \sF$, which leads to the troublesome relation in Eq. (P). This is a general result 
irrespective of whether fi is generated spontaneously or radiatively. One cannot get both \x 
and J B^ of the order of Mz and EWSB fails because the large value of B^ destabilizes the 
scalar potential. 

Motivated by the NMSSM model, we consider a minimal gauge-mediated model with an 
extra singlet Higgs field. The most general superpotential respecting a Z 3 symmetry is 

W = £ S S$$ + £ N N$$ - ^-S 2 N - ^-N 2 S + X s SH u H d + \ N NH u H d - ^N 3 . (7) 



The gauginos and scalars acquire mass radiatively at the 1-loop and 2-loop level respectively, 
and have the following form 

cv ■ 

M t ~ h-^A, (8) 

47T 

m)^2J:c{k l (^A)\ (9) 

where k x = 5/3, k 2 = k 3 = 1, and Cf = V 2 for U(l), c{ = 3/4 for weak SU(2) doublets 

and Cl = 4/3 for color triplets. A = F$/S is the effective SUSY-breaking scale. In models 
in which messengers and matter do not interact directly, the trilinear soft A-terms and 
m 2 N arise only at 2-loop and 3-loop respectively, and can be neglected. The superpotential 
Eq. (0) contains messenger-matter couplings, thus inducing A\,k and m 2 N at 1-loop and 
2-loop respectively M, 



m 



r^j 



7 8 3 13 \ / A \ 



10« €jv ( -a (N - ~a 3 - ~a 2 - —a ± j ^— J , (11) 



where a^ N = t-^/Aix. Note that the contributions from the gauge couplings can make m% 
negative for small values of £jv- Among the couplings in Eq. ([?]), we may anticipate that 
Cs, £n, Xn, k ~ Oil) and that t/n, rjs and As will be small due to the large VEV from S. We 
will determine the pattern for phenomenologically viable solutions. The only result contrary 
to these expectations is that we find r] N ~ 0(1). We will not explore the possible symmetries 
in terms of model-building that naturally lead to this pattern |6|J^,|lO[1 . 
The relevant part of the scalar potential is 







Vlfiggs — Vp + Vo + Vtadpole + Vsoft + AV, 




(12) 


where 






V F = 


X S S + X N N\ 2 (\H U \ 2 + \H d \ 2 )+ \X N H u H d - kN 2 -^-S 2 - Vn NS\ 2 




+ 


IXsH^ + Fs-^^-rjsNS] 2 , 


(13) 


v D = 


9 l(HiaH u + H d aH d ) 2 + 9 ^(\H u \ 2 -\H d \ 2 ) 2 , 


(14) 


^tadpole 


ZsZnNF 2 
8tt 2 S ' 


(15) 


V so ft = 


= m 2 Hu \H u \ 2 +m 2 Hd \H d \ 2 +m 2 N \N\ 2 -(X N A XN NH u H d + ^N 3 + h.c.), 


(16) 


AV = 


3 

327T 2 


f^r 3 m 2 3 m 2 3 
mUln * )+mj (In I 2 ) 2m*(ln * ) 

tlV Q 2 2 J t2K Q 2 2 ! tK Q 2 2 J 


> 


(17) 



where Vt adpo ie arises from the N<&<& coupling via one-loop tadpole diagrams involving the 
messenger fields [OT. For the one- loop potential term AV, we have only included the dom- 
inant contributions from the top-quark and top-squark diagrams. Q is the renormalization 
scale in the DR scheme. 



III. ANALYSIS OF THE MODEL 

A. Numerical Procedure 

To study physics at the electroweak scale, we perform a one-loop renormalization group 
(RG) evolution analysis. The relevant RGEs can be derived from fll] and are listed in 
Appendix A. 

We fix the boundary conditions for the RG evolution as follows: At the messenger scale, 
we impose the conditions in Eqs.(8)-(11). The coupling constant £s is fixed at the messenger 
scale to be unity. £,n,Vs an d t]n are solved at the electroweak scale by minimizing the 
potential as described below. Values of the other three couplings, As, A^ and k are chosen 
randomly at the electroweak scale. We require the couplings to be in the perturbative ranges 

- 1 < £n(S) < 1, < A s , Ajv, »7jv < 1, < k < 0.65, (18) 

where the upper bound on k comes from the demand that it not hit the Landau pole below 
the GUT scale. If we require that k be perturbative only below the messenger scale A where 
unspecified new physics sets in, we find k < 1.36. It is found to be unnecessary to place 
bounds on the value of rjs- In the numerical analysis, we will neglect Yukawa coupling terms 
in the RGEs, except for those involving the top quark (h t ). The A-terms run slowly and so 
can be approximated by their boundary values. All RGEs are run down to Q = ^/m^m^. 

In the RGE running, m 2 Hu receives large correction from the heavy stop loop through the 
Yukawa coupling h t , 

m 2 Hu (Q) ~ m 2 Hu (S) - ^h 2 (ml(S) + m 2 iR (S) + m 2 Hu (S)) In- (19) 

and becomes negative at the electroweak scale, triggering EWSB. A satisfactory estimate 
for m 2 N at the electroweak scale can be found by approximating the RGE for m 2 N : 



m 



%(Q) ~ m%(S)-^- 2 (2m 2 Hd {S) - -^- 2 h 2 t m}{S) In | + m%(S) + A\ 




[2 (3m%(S)+Al)ln^. (20) 

which can be either positive or negative depending on the relative values of the couplings. 
The physical mass spectrum at the electroweak scale is determined by the five VEVs 

(H u )=v u , (H d )=v d , (N)=x, (S}=y and F s . (21) 

The values of Mz and top-quark mass fix the scale 



v = Jv 2 u + v 2 d = 174 GeV, m t = h t v u = 165 ± 5 GeV (22) 



where we have used the DR value of the top-quark running mass. We also require that 
tan/? = v u /vd satisfy 

1.5 < tan/3 < 50, (23) 



which leads to 0.92 < h t < 1.2. We impose the following constraints on the other VEVs 

100 < x < 1000 GeV, 3.8 x 10 4 < A, y < 10 5 GeV, A < 0.9y, (24) 

where A = Fs/y. The lower bound on A arises from the right-handed slepton mass constraint 
(see below), and the last condition in Eq. (|24]) is adopted to avoid possible large corrections 
from A p« y. The calculated mass spectrum for the sparticles must satisfy the current lower 



limits from direct experimental searches ||T2"|.|I"3" 



nig > 190 GeV, m~ eR > 80 GeV, 



m h = m Sl > 95.5 GeV, m^ = mpj > 84.5 GeV, m H ± > 77.3 GeV, 
m-± > 95 GeV, m- h > 75 GeV, 



m^,o m.oo 



if _*i < i m > 55 QeV; if — ^ > 1, m fl > 68 GeV; 

if — ^ > 1 and m^ < 87 GeV then m fl > 84 GeV. (25) 

The minimization conditions to be imposed at OiMz) can be derived by differentiating 
Eq. flTJD with respect to v d , v u and x: 

—iii i\ 
2m 2 Hd v d H — | — + w u Aiv {-2kx 2 - 2xyr] N - y 2 r] S + 2v d v u X N ^ + 2w d (xAAr + y\ s f 

+v u X s (2F S - x 2 VN - 2xyi] S + 2v d v u X s ) + d Vd AV = 0, (26) 

Til) (ii ?j i 

2m\v u ^— | — + v d \ N \-2kx 2 - 2xyrj N - y 2 r] S + 2v d v u \ N ) + 2v u {x\ N + yX s f 

+v d \ s (2F S - x 2 r] N - 2xyr] S + 2v d v u X s ) + d Vu AV = 0, (27) 

2m 2 N x + s N S + (2kx + yrj N ) hkx 2 + 2xyrj N + y\s - 2v d v u X N ) + 2v 2 X N (x\ N + yX s ) 
87v z y v ' 

+ (xr] N + yrjs) (-2F S + x 2 7] N + 2xyr] S - 2v d v u X s ) + d x AV = 0, (28) 



where g 2 = g' + g 2 . Running £s down to the electroweak scale, we get £,s(Q) — 0.94. r]s, tjn 
and £n are determined from the minimization conditions Eqs. (p6|)-(p8|). These equations 
also lead to the relations 



M\ , [m 2 Hd + (d Vd AV)/(2v d ) ] - [m 2 Hu + (d Vu AV)/(2v u )] tan 2 (3 

tan 2 (3 - 1 



^ = _^Z + l-^ v^d-w/v-^j v-n u ,-„„-.„ V ~„, J ..~~ ^ (29) 



Sm W = [m 2 Hd + (d Vd AV)/(2v d )] + [J Hu + (d Vu AV)/(2v u )} + 2/x*' (30) 

where 

(i = X N x + X s y, (31) 

2 2 
B^ = X n (A Xn x + kx 2 + r] N xy + -^ X N v u v d ) + X s (-F s + r] S xy + -^ A^u,*). (32) 

The parameters /i and B^ can be evaluated by either set of equations and they must be 
consistent with each other since all the parameters have been determined by conditions on 
EWSB. 
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B. Viable solutions, /j, and B, L parameters 
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FIG. 1. Phenomenologically viable solutions in the \s~£ N , Xs-rj s , Xs~r] N an( ^ A,g-tan/? planes. 

On inspection of the parameter space of the variables chosen randomly in the allowed 
ranges in Eq. flI8|) , it is found that there are solutions satisfying the theoretical and exper- 
imental constraints of Eqs. (p4D-(p5|). We first note that rj^ and the "NMSSM couplings" 
Xn and k are essentially independent of other parameters and mostly of 0(1), 



0.6 < 7] N < 1.0, 0.05 < X N < 1.0, 0.001 < k < 0.65. 



(33) 



£n and r/s depend quite crucially on the values of Xs and take on values of 0(1) and 
C(10" 3 - 10~ 2 ) respectively, 



corresponding to 



0.95 < ^ < 0.95 and |^|>0.1, (34) 

0.01 < r]s < 0.05 and \r] S \>W-\ (35) 

< A s < 4 x 1(T 3 . (36) 



The correlations for those parameters are shown in Fig. [|. There are two rather distinc- 
tive regions in \g. 

1. Vanishing As: There are viable solutions for vanishing \$. In this case, typically 
As ~ 10~ 3 , leading to tightly constrained solutions for the allowed values of £at and rjs- 

- 0.95 <£ N < -0.3, -0.01 < T] S < -4 x 1(T 3 . (37) 

It is interesting to note that rjs and £jv are always negative, thus picking out a particular 
choice of phases in the superpotential 

W D £ S S$$- \£ N \ A$$ + ^-S 2 N + ... 

Also notice that small tan/5 is excluded. 

2. Finite As: For As ~ 10~ 3 , the couplings ^ and r\s can take on the whole range of 
permissible values. Moreover, small tan/5 solutions are allowed. 

Since we wish to keep all terms in the superpotential that respect the Z3 symmetry, we 
avoid a vanishing As and restrict our attention to the case As > 10~ 3 in the rest of the 
paper. This choice gives us greater freedom in parameter space. 

The parameters \i and B^ are spontaneously generated and are of 0(M Z ) and (9(M|), 
respectively. As shown in Fig. 0(a), the values lie in the ranges 

275 < fi < 550 GeV, 5 x 10 3 < B^ < 2 x 10 5 (GeV) 2 . (38) 

Moreover the "5-parameter" (B^/n) familiar from supergravity theories is also 0{Mz) for 
most of the parameter space 

10 < -± < 300 GeV, (39) 

as shown in Fig. |2|(b). For tan/5 > 2, M|/2 ~ — fi 2 — m 2 H . A good measure of fine-tuning 
of the Higgs potential parameters is Mf /2/z 2 ||14j| . The fine-tune is at the 1-6% level as seen 
in Fig. |3], where the fine-tune parameter is shown versus tan/5 and A. This level of fine-tune 
could be significantly improved if one introduces extra vector-like triplets that couple only 
to the singlet N |TjJ. Recall that in the MSSM, fi of electroweak scale is put in by hand. 



Even if the very lowest value of \x (consistent with mg R > 80 GeV) is used, the best possible 
fine-tune is 7%. 

8 



> 

CD 
O 



CQ 



10- 



10" 





1 ' ' ' ' 1 ' ' ' 


i 1 i i i. 


- 




- 


- 


,">'■'--. ' 


- 


" 


'•*'"' •' 


■ 




:^ '::,V- 








- 


■/ V''- --" : • ■■"., - 


- 


" 


!,,,,!,,, 


, 1 , , ," 



300 



400 
fj, (GeV) 



500 



> 

CD 
O 



m 





1 




1111 1 


1 1 


200 




% 




— 


100 






— 










- 








T©^ 




50 






'"%;: ^ 


- 








•■»V-«---.--. 


- 








■■'.:■'■■ ■■ ■ 


- 


20 


1 




1 , , , , 1 1 


. 1 



10 
tan/S 



20 



50 



FIG. 2. The allowed region in the (a) B^-fi plane and (b) (B M //i)-tan (3 plane. 



o 
o 



a. 



CD 

a 

CD 

a 




50 






60 



70 
A (TeV) 



80 



90 



FIG. 3. The percentage fine-tune versus tan/3 and A. 



C. Mass Spectrum 

The details of sparticle spectroscopy of models with gauge mediation have been discussed 
in the literature over quite large ranges of both tan (3 and A [|1^]. Many of those results carry- 
over unchanged to the case at hand. As one would expect, the main differences arise from 
the specific choice of the Higgs sector superpotential. The mass-squared matrices for the 
CP-even and CP-odd Higgs bosons and their eigenvalues are presented in the Appendix B. 

We find the range for the CP-even Higgs boson masses 



95.5 < m Sl < 130 GeV, 200 < m s , 2 < 650 GeV, m S3 > 5 TeV. 



(40) 



We see that m Sl has the usual upper bound of 130 GeV. If the perturbativity for k is relaxed 
to the messenger scale, m,s 1 can be as heavy as about 170 GeV. In the NMSSM, there are 
two CP-odd Higgs boson mass eigenstates (Pi,P 2 ). The mass of the lighter CP-odd state 
lies in the range 340 < m Pl < 650 GeV indicating that we are well into the decoupling 
regime Jl7|]. Pi is primarily a mixture of ImP^ and lmH u since the mixing angle 7 between 
A and ImN is O(10~ 4 ) (see Appendix B). We also find that the charged Higgs bosons H ± 



are almost degenerate with Pi. 

Since /i is large compared to Mz, the higgsinos and gauginos decouple from each other 
and the lightest neutralino is mostly bino. It is the x?-NLSP for a large region of parameter 
space 



tan/? < 22 or A < 64 TeV. 



(41) 



Due to the mixing in the r-slepton (stau) mass matrix, m fl can be much smaller than mg Rt jx R 
and becomes the NLSP for large tan/3 and high values of A. These features are shown in 
Fig. f|. We find the mass ranges for the lighter sparticles are 

65 < rn^o < 125 GeV, 120 < m^+ < 230 GeV, 68 < m fl < 150 GeV. (42) 

Gluinos and squarks are much heavier than the electroweak states because they acquire 
mass mainly from strong interactions. We find the ranges the of stop masses as follows: 

650 < m h < 1080 GeV, 700 < m h < 1140 GeV. (43) 

Table | gives a representative set of parameters and the corresponding mass spectrum. 



tan/3 


£s 


6v 


Vs 


VN 


As 


Xn 


k 


(S) 


A 


(N) 


24.3 


0.94 


0.71 


8.3 x 10" 3 


0.87 


3.0 x 10" 3 


0.31 


0.40 


97 TeV 


77 TeV 


413 




1* 


Mi 


M 2 


M 3 


m x? 


m xt 


msi 


m Pl 


m H ± 


m fl 


17lf 2 


m h 


m h 


m o 


421 


107 


208 


649 


106 


196 


122 


500 


505 


109 


288 


949 


982 


10.7 TeV 



TABLE I. A representative set of parameters and the mass spectrum. All masses are in GeV 
unless specified. 
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IV. MESSENGER SNEUTRINOS AS COLD DARK MATTER 

The possibility of messenger sneutrinos (0°) as a cold dark matter (CDM) candidate was 
investigated in the context of the minimal GMSB model in Refs. |l8|Jl^1 . Requiring that the 



Universe not be over-closed by the messenger sneutrinos (fi^oh 2 < 1) puts an upper bound 



on their masses 



m,j.o < 3 TeV. 



(44) 



It certainly requires some degree of fine-tuning to obtain such a light messenger particle as 
the messenger scale is about 50 — 100 TeV. More severely, a scalar DM particle with mass 
lighter than 3 TeV and with SM interactions has already been experimentally ruled out ||20|| . 
In the GMSB model with an extra singlet Higgs field, additional contributions from the 
singlet N can sufficiently increase the annihilation cross-section of 0°, and therefore loosen 
the above bound. 

We perform a relic density calculation in our model. We first note that the freeze- 
out temperature Tj ~ m^o/20 ^> m^- — m^o, which implies that both charged and neutral 
components are present in the thermal bath at Tf | ID| . We calculate the relic density in the 
limit of symmetric phase for the four components ((ffi , <fi^ with Q = 0, — ). On examining 
the interaction vertices, by far the dominant annihilation processes are 



NN, H„H, 



(45) 



via the t- and w-channel messenger exchange and s-channel singlet exchange. Following ||21|| . 
we define a Lorentz invariant function 



w(s) 



^rtlll |.M(<W - ab)\Hcos9 

Q,Q,a,b 



(46) 
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FIG. 5. Messenger sneutrino masses that do not over-close the universe, 
where \fs is the center of mass energy. Specifically in our model, 



w(s) 



f( VN k(S)) 2 $ 



2 r 



y/2 ^ k ' 3 I s / 



where d = 2£jv(2£s — r] N )/(r] N k). With k = Tf/m^o and 



tf 2 



Ami 



A = w(4mj ), B = 3 



■2ml -^ - w{s) 



ds 



s=Am? r 



the thermal average (cn> re i) can be expressed as 

1 



(<™rel) 



mi 



(A + Bk) 



(47) 



(48) 



(49) 



It is customary to express the relic abundance in terms of the mass density in units of the 
critical density Q^o = p^j ' p c . It is found that 



a 



2 _ 8.5 x 1(T 5 /m o\ 



K 



VTeV/ A+±B 



1d k ' 



(50) 



where g* ~ 232.5 for the particle content of the NMSSM at the freeze-out temperature Tf. 
The value of k is related to Tf and is obtained iteratively by 



k = m^o /Tf = In 



0.076 M Pl 



(A + Bk)^[k 



(51) 
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From the parameter space that provided a phenomenologically viable solution to the \x- 
problem we select the subspace that does not over-close the Universe. Fig. |5] shows the relic 
density versus the cold dark matter mass m^o, where m^o/(S) was allowed to vary between 
0.1 and 1. If we assume that Q^oh 2 ~ 0.3, then the mass of the lightest messenger sneutrino 
lies in the range 

6 < m o < 25 TeV. (52) 

The smaller values of m^o correspond to m^o/(S) close to 0.1 and are undesirable from 
the fine-tuning argument. We conclude that massive scalar messenger cold dark matter is 
natural in this model. 



V. A RESOLUTION OF THE COSMOLOGICAL DOMAIN WALL PROBLEM 

It is well-known that the NMSSM predicts the formation of domain walls at the weak 
scale arising from the spontaneous breaking of the Z 3 symmetry of the superpotential. Al- 



though there are proposals to resolve this problem ||22|| , the most straightforward solution 
still involves the introduction of higher dimensional non-renormalizable operators that ex- 
plicitly break the discrete symmetry. It has been shown that this solution is unsatisfactory 
when SUSY-breaking is mediated by gravity [|J . We briefly indicate where the problem lies. 

A wall network propagates once it has formed. The most important force acting on the 
domain walls that can dissolve them is pressure. Pressure arises from an explicit violation 
of the discrete symmetry. If the vacua are separated by an energy density e, the pressure is 
typically of order e. When this difference is greater than the surface tension a/R, where R is 
the curvature scale of the walls and a is the surface energy density, there is a minimum value 
of e above which the dynamics will have been dominated by pressure before today. Thus, 
as the Universe expands, pressure catches up with surface tension, ultimately dominates it, 
the wall network is destroyed and the true vacuum is selected. 

In general if the Z3 symmetry is broken at an energy scale vo, then a ~ vq 3 , and if 
pressure is to be dominant at an epoch T* before the quark-hadron phase transition, the 
curvature scale will be roughly the scale factor R ~ Mpi/T 2 , and the relation 

6 ~ i^r (53) 

must hold. With a real scalar field having a VEV Vq, the leading higher dimensional 
Z 3 -breaking operator can be written as \'<j) 5 /Mpi. Assuming other couplings to be 0(1), 
the above constraint can now be rephrased as a constraint on A': 

A<>(!) 2 . (54) 

Note that the above relation holds only when all the fields in the Z3-breaking operator 
acquire the same VEV as the scale at which the symmetry is broken. The relevance of this 
point will become clear below. 

Nucleosynthesis dictates that the walls decay before the Universe cools to 0.1 to 1 MeV. If 
this were not the case, the abundances of the lighter elements would be significantly smaller 
than measured. With T* ~ 10 MeV and v = v ~ 174 GeV, 
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A' > 1(T 8 . (55) 

However, when the mediation of SUSY-breaking is gravitational, this leads to destabilization 
of the hierarchy [f§||. For example, the operator \' (H u H d )N 2 / M Pi along with X N H u H d N 
generates a tadpole that leads to a term in the superpotential of the form fi'N with 



P' ~ nr„-?\9 A Gravit V i ( 56 ) 



where A Gravity ~ O(10 w GeV) is the SUSY-breaking insertion. For A' > 1CT 8 , this gives 
(/i'x) 1 / 3 > 6 TeV thus destabilizing the electroweak scale. In other words, to obtain 
fjfx ~ 0(v 3 ), one needs A' ^ 10~ n , in obvious contradiction with cosmology. 

In our model, the operator that can destabilize the hierarchy most severely is 
A' \H u Hd)S 2 1 'Mpi. The only field with a messenger scale VEV is S. The lower bound 
on A' is now altered to 

A'A 5 >(^) 2 ~HT 3 (57) 

for v = y ~ 10 5 . With A s ~ 0(l(r 3 ), this leads A' to the natural value A' ~ 0(1). 
Such a large value of A' looks a lot worse than in the gravity-mediated case in terms of the 
destabilization of the hierarchy. However, a crucial observation is that in GMSB models the 
SUSY-breaking scale is typically very low. In the model under consideration, SUSY-breaking 
is communicated via a direct interaction £ S SQQ and consequently the fundamental scale of 
SUSY-breaking F = F$/£,s ~ A 2 ||, where A is 50-100 TeV. As earlier, the operators 
A' \H u Hd)S 2 1 'Mpi and XsH u HdS generate a term in the superpotential of the form fi'S. In 
this case, 

A*' ~ T^WA 2 ~ f-^O (—^—) 2 (30 GeV) 2 . (58) 

p (16vr 2 ) 2 \4x 10- 3 ; VIOOTeV/ K ' K ' 

Therefore, the stability of the electroweak scale can be maintained, fi'y ~ 0(v 3 ), while 
simultaneously satisfying the cosmological constraints. However, if the fundamental SUSY- 
breaking scale happens to be much higher Fq 3> A 2 , the problem of destabilizing the elec- 
troweak scale will remain. 



VI. DISCUSSION AND CONCLUSION 

By design in our model, both \i and B^ parameters are generated spontaneously at the 
right order, thus solving the /z-problem. The fine-tune for the cancellation between /z 2 and 
m 2 H to obtain the correct Mz is still at the percentage level, as in the MSSM. The primary 
reason that we are able to overcome the /i-problem in GMSB || is the inclusion of messenger- 
matter interactions in the superpotential. However, the smallness of the parameters rjs 
and As should be further explored in terms of discrete symmetries arising from the SUSY- 
breaking sector. 
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As for low-energy phenomenology, we see from Table p] that the sparticle spectrum is 
typical of gauge-mediated models. It is worth mentioning that if we require k to be pertur- 
bative only below the messenger scale A, we find k < 1.36. This results in the upper bound 
on the lightest Higgs boson mass to be about 170 GeV. 

In summary, we studied in some detail the next-to-minimal super symmetric standard 
model with the gauge mediation of supersymmetry breaking. We found that it is feasible to 
spontaneously generate the /i and B^ terms which are consistent with electroweak symmetry 
breaking and the current sparticle mass limits. Messenger sneutrinos with mass in the range 
6 to 25 TeV can serve as cold dark matter. It is also possible to find a solution to the 
cosmological domain wall problem consistent with EWSB due to the low scale of SUSY- 
breaking 
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APPENDIX A: 

The 1-loop renormalization group equations of the couplings and soft-SUSY breaking 
terms relevant for our analysis are listed. We only include contributions from the third 
generation to the superpotential and soft-terms. 

W = £ S S$® + 6vA$<5 - ^-N 2 S - ^-S 2 N + X N NH u H d 

+ X s SH u H d - ^N 3 + h t H u QU + h b H d QD + h T H d LE } (Al) 

V soit = m 2 H \H d \ 2 + m 2 H jH u \ 2 + m 2 N \N\ 2 + m 2 j£ L | 2 + m 2 iR \t R \ 2 + mljb R \ 2 + m 2 Jf L | 2 

+ m 2 R \f R \ 2 + -M X BB + -M 2 WW + -M 3 ~gg + {A t h t t L t R H u + A b h b b L b R H d 

k 
+ A T h T f L f R H d - X N A XN H d H u N - -A k N 3 + h.c). (A2) 



For the gauge and Yukawa couplings: 

lGvr 2 ^' = 1 V 3 , l^ 2 j t 92 = gl 16tt 2 ^ 3 = -3^, (A3) 

l^ 2 jht = h t (6h 2 + h 2 b + X 2 N + X 2 s - fg' 2 - 3g 2 2 - ^) , (A4) 

l&r 2 ^ = h fa + h 2 + h 2 r + X 2 N + X 2 s - \g' 2 - 3gj - fgl) , (A5) 

16tt 2 — h T = h T (Ah 2 T + 3h 2 b + X 2 S - 3g' 2 - 3gfj . (A6) 
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For the couplings particular to the NMSSM and to our model: 

l^ 2 jk = 6k(k 2 + \ 2 N )-r ] 2 N (2k + r ]s ), (A7) 
lQ7T2 Jt XN = Xn ( 2A ^ + r >s + r >N + 5 d + ^n + 2A; 2 + Sh 2 + 3h 2 b +h 2 T - g' 2 - 3g 2 2 ) , (A8) 

16n2 Jt Xs = Xs ( 2A ^ + 4A ^ + r >s + ri 2 N + ^s + 3h 2 t + 3h 2 b + h 2 T - g' 2 - 3«? 2 2 ) , (A9) 

16tt 2 ^6v = Zn (9& + 4£ 2 + 77I + V 2 n + 2A 2 V + 2k 2 + ^-g' 2 - 3^ 2 - y^ 2 ) , (A10) 

16tt 2 ^ = £ S (9^ + 4& + „* + < + 2A| + f</ 2 - 3^ 2 - ^g 2 ) , (All) 

16tt 2 ^ W = ^ (5^| + 10& + 3 V 2 + 2,7)1 + 2A| + A\ 2 N + 4k 2 ) , (A12) 

1 ^ 2 j t Vs = Vs (5^ + ml + Zvl + 3< + 4A 2 S + 2X 2 N + 2k 2 ) . (A13) 

For the trilinear couplings of the third generation: 

-2 « 4 _ iol,2 4 OI.2J , o\2 /l /i ^^ /2 3 2 8 2 



IGtt^-A = 12^A + 2/z 2 A 6 + 2A^ A - 4 (^"Mi + -<fe 2 M 2 + -g 2 3 M 3 ) , (A14) 

16tt 2 ^A 6 = l2h 2 b A b + 2h 2 A t + 2h 2 T A T + 2A^ A - 4 {^g' 2 M l + \g\M 2 + j^f M 3 ) , (A15) 

16tt 2 -A t = 8/i 2 A r + 6h 2 A b + 2A 2 V A A - 6 (g' 2 M 1 + ^ 2 2 M 2 ) , (A16) 

16tt 2 ^A a = 8A^ A + 2A|A A - Ak 2 A k + 6/i 2 A t + 6h 2 b A b + 2h 2 T -2 (g' 2 M x + 3g 2 2 M 2 ) , (A17) 

167r 2 ^A fc = 12 {k 2 A k - A^ A ) + r, 2 N A k . (A18) 
For the soft masses: 

l6n 2 -M 1 = 22g f2 M 1 , \§it 2 jM 2 = 2g 2 2 M 2 , 16tt 2 ^M 3 = -6^ 3 2 M 3 , (A19) 
16vr 2 -m 2 L = 2/i 2 (m? L + m\ u + m\ R + A 2 ) + 2h\ (m\ + m\ d + m\ R + A 2 

~/2/- o /" „/2 71 t2 1 ^ 2 2 1 \2./2 



+ ^g'X - 8 (ggr^i + 402M2 + 3^ 3 M 2 J , (A20) 

l^ 2 j t ml = Ah 2 (ml + m 2 Hu + m\ R + A 2 ) - 8 (^' 2 M 2 + ^ 3 2 M 2 ) - ^' 2 C , (A21) 

167r 2 |m? fl = Ah 2 b (m| 3 + m| d + m? + A 2 ) - 8 (±g' 2 M 2 + ^ 2 M 3 2 ) + ^' 2 (, (A22) 

16vr 2 |m 2 L = 2/> 2 « + < + < + A 2 ) - 8 Q,' 2 M 2 + ^ 2 M 2 ) - </ 2 C, (A23) 

^n 2 jm 2 R = 4/i 2 (m 2 L + m\ d + m\ R + A 2 ) - 8g' 2 M 2 + 2g' 2 (, (A24) 
lfa 2 -m 2 Hd = 6h 2 b (m 2 Hd + m 2 iL + m\ R + A 2 ) + 2A| (m 2 ^ + m^) 
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+ 2/4 { m H d + ™% L + m\ R + 4) + 2X 2 N (m 2 Hd + m 2 ^ + m 2 N + A 



~ 8 (\ 9 ' 2 M 2 + ^ 2 2 M 2 2 ) - </ 2 C, (A25) 

-8 Q«/ 2 M 2 + ^ 2 M 2 ) + 2X 2 N (m| d + m| u + m 2 , + A 2 ) + </ 2 C, (A26) 
levr^m 2 , = 4A 2 , (m 2 , + m 2 Hu + m| d + A 2 ) + 4k 2 (3m 2 , + A 2 ) , (A27) 

where ( = Y^i Yim 2 is the hypercharge-weighted sum of the squares of the soft masses. As a 
consequence of the boundary conditions at the messenger scale, ( = at all energies. 



APPENDIX B: 

The squared mass matrices for the neutral and charged Higgses are presented. The 
analytic expressions for their diagonalization are displayed. 

where v\ = v d , v 2 = v u , v 3 = x and 

a = — (g 2 vl + v u ((2kx 2 + 2xA Xn + 2xyi] N + y 2 i] S )X N + (~2F S + x 2 r] N + 2xyr] S )\s)) 

-(d Vd AV)/v d + d 2 d>Vd AV (B2) 

b = —(fvl + v d ((2kx 2 + 2xA Xn + 2xyr) N + y 2 r] S )\ N + (-2F S + x 2 t] N + 2xyr) S )X s )) 







-(ft.A^K + e.AV (B3) 



., v -16£;7rVy,4 fe - F 2 ^ s + IGtt^ 2 ^^ + 8ir 2 y(y VN (6kx 2 + (3x 2 - y 2 ) Vs 
8n z xy 

+2v d v u \ N ) + 2x^r] 2 N + 2(4k 2 x 3 + A Ajv v d f u A^ + yv d v u r) S \s - yvjX N X s - yvl\ N \ s ))) 

-(d x AV)/x + dl x AV (B4) 

d = -{2kx 2 + 2xA Xn + 2^77^ + y 2 r] S )\ N + A 5 (2F 5 - rr 2 ?^ - 2xyr) S + 4^ u As) 

-^UdUu + ^^A^ + ^^^Ay (B5) 

e = 4w d AAr(a;AAr + yAs) - 2v u ((2kx + A Ajv + yr) N )\ N + (xr) N + yr) S )\ s ) + d 2 VdX AV (B6) 

/ = 4v u X N {xX N + y\ s ) - 2v d {{2kx + A Xn + yrj N )X N + {x Vn + yvs)X s ) + d 2 uyX AV. (B7) 

The eigenvalue equation from this matrix is 

(m 2 f + a 2 (m 2 ) 2 + ai m 2 + a = (B8) 

where 
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With 



a = -abc - 2def + af + be 2 + cd 2 (B9) 

a 1=a b + ac + bc-{d 2 + e 2 + f 2 ) (BIO) 

a 2 = -(a + b + c). (Bll) 



1 1 1 1 1 r 

-(a 1 a 2 -3a ) - — a\, q = --oi + -o 2 , 6> = - arccos ^r^, (B12) 



the eigenvalues are 



m 



?7i 



^(cos + V3 sin 0) - — (B13) 

o 

^(-costf + v^sintf)- — (B14) 



m| 3 =2^cos^-y. (B15) 

Now for the CP-odd Higgs bosons. The mass squared matrix is 

M 2 CP . odd = (J J s ) (B16) 

where 
r = 2^1 (( 2A;a;2 + 2xA a jv + 2 ^^ + 2/V) Ajv + (~2F S + x (xr^ + 2y Vs )) \ s ) 

s = w^l 4Skn2x2 y A k - FsZnSs + 167r 2 yF 5 (2x77^ + VVs) - 8n 2 y{r] N {2kx 2 y + y (x 2 + y 2 ) Vs 
-2v d v u (y\ N + 2xA 5 )) + 2 (yrjs (2kxy - v d v u \ s ) + \ N ((-4&x - A Xn ) v d v u + yv 2 X s )))} 

j = -(y ((2kx - A Xn + yr) N ) \ N + (xrj N + yrj s ) X s )) . 

The transformation matrix that takes us from the (A = sin/^Imif^ + cos / 3Im.£/ u ) — ImiV 
basis to the diagonal basis is 



cos 7 sin 7 
- sin 7 cos 7 



(B17) 



where 



sin2 7 = . 2 l cos2 7 = . S L (B18) 



and the eigenvalues are 



s — r) + 4j 2 J (s — r) + 4j 



"4,*. = o ( s + r =F V ( s " r ) 2 + 4 ^)- ( B19 ) 



Finally, the squared mass of the charged Higgs is 

m 2 H± =M$y- (X 2 N + X 2 s )v 2 + r (B20) 

where r is the first entry of the CP-odd Higgs matrix. 
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